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Abstract 



d 

We derive a forward partial integro-differential equation for prices of 
call options in a model where the dynamics of the underlying asset under 
^ , the pricing measure is described by a -possibly discontinuous- semimartin- 

' gale. This result generalizes Dupire's forward equation to a large class of 

QQ I non-Markovian models with jumps. 
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Since the seminal work of Black, Scholes and Merton [51 US] partial differen- 
tial equations (PDE) have been used as a way of characterizing and efficiently 
computing option prices. In the Black-Scholes-Merton model and various exten- 
sions of this model which retain the Markov property of the risk factors, option 
prices can be characterized in terms of solutions to a backward PDE, whose 
variables are time (to maturity) and the value of the underlying asset. T he use 
of backward PDEs for option pricing has been extended to cover options with 
path-dependent and early exercise features, as well as to multifactor models 
(see e.g. [T]). When the underlying asset exhibit jumps, option prices can be 
computed by solving an analogous partial integro-differential equation (PIDE) 

A second important step was taken by Dupire [16l [171 US] (see also Dcrman 
& Kani [15 J who showed that when the underlying asset is assumed to follow 
a diffusion process 

dSt = a{t, St)dWt 

prices of call options (at a given date tg) solve a forward PDE in the strike and 
maturity variables: 

dT ^ ' dK 2 dK-^ 

on [to, oo[x]0, oo[ with the initial condition: V/C > Ct^ {to, K) = {Sto — K)-^-. 
This forward equation allows to price call options with various strikes and ma- 
turities on the same underlying asset, by solving a single partial differential 
equation. Dupire's forward equation also provides useful insights into the in- 
verse problem of calibrating diffusion models to observed call and put option 
prices 

Given the theoretical and computational usefulness of the forward equation, 
there have been various attempts to extend Dupire's forward equation to other 
types of options and processes, most notably to Markov processes with jumps 
[H [ini [121 [23 [3] ■ Most of these constructions use the Markov property of the 
underlying process in a crucial way. 

As already noted by Dupire [18], the forward PDE holds in a more general 
context than the backward PDE: even if the (risk-neutral) dynamics of the 
underlying asset is not necessarily Markovian, but described by a continuous 
Brownian martingale 

dSt = StatdWt 

then call options still verify a forward PDE where the diffusion coefficient is 
given by the local (or effective) volatility function a{t, S) given by 

a{t,S) = ^E[a^\St = S] 

This method, also known as "Markovian projection" , is linked to the construc- 
tion of a Markov process which mimicks the marginal distributions of a martin- 
gale [1[71[51[17]. 
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We show in this work that the forward equation for call prices holds in a 
more general setting, where the dynamics of the underlying asset is described by 
a -possibly discontinuous- semimartingale. The derivation of this result is based 
on the Tanaka-Meyer formula and has the merit of not requiring ellipticity or 
non-degeneracy of the diffusion coefficient; the results are thus applicable to pure 
jump models and point process models used in equity and credit risk modeling. 

Our result extends the forward equation from the original diffusion setting 
of Dupire jl7j to various examples of non-Markovian and/or discontinuous pro- 
cesses and implies previous derivations of forward equations [H [TUl [51 HH [HI [TSl 
[23 [55] as special cases. Section [5] gives examples of forward PIDEs obtained 
in various settings: time-changed Levy processes, local Levy models and point 
processes used in portfolio default risk modeling. In the case where the under- 
lying risk factor follows, an Ito process or a Markovian jump-diffusion driven by 
a Levy process, we retrieve previously known forms of the forward equation. In 
this case, our approach gives a rigorous derivation of these results under precise 
assumptions in a unified framework. In some cases, such as index options (Sec. 
12. 5p or CDO expected tranche notionals (Sec. 12. 6p . our method leads to a new, 
more general form of the forward equation valid for a larger class of models than 
previously studied [31 [TH [51] . 

The forward equation for call options is a PIDE in one (spatial) dimension, 
regardless of the number of factor driving the underlying asset. It may thus 
be used as a method for reducing the dimension of the problem. The case of 
index options (Section 12. 5|) in a multivariate jump-diffusion model illustrates 
how the forward equation projects a high dimensional pricing problem into a 
one-dimensional state equation. 

1 Forward PIDEs for call options 

1.1 General formulation of the forward equation 

Consider a (strictly positive) price process S whose dynamics under the pricing 
measure P is given by a stochastic volatility model with jumps: 

/•T /'T /'T /'+00 

St^So+ r{t)St-dt+ St-dtdWt+ / St-ie^ - l)M{dt dy) {!) 

Jo Jo Jo J-oo 

where r{t) > represents a (deterministic) bounded discount rate, 5t the (ran- 
dom) volatility process and M is an integer- valued random measure with com- 
pensator fj,(dtdy;Lu) = 171(1, dy, Lu) dt, representing jumps in the log-price, and 
M = M — /z is the compensated random measure associated to M. (see [13] for 
further background). Both the volatility 6t and m{t,dy), which represents the 
intensity of jumps of size y at time t are allowed to be stochastic. In particular, 
we do not assume the jumps to be driven by a Levy process or a process with 
independent increments. 

We assume the following integrability condition: 
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Assumption 1. 



VT > 0, E 



2 Jo 



exp 



< oo (H) 



The value Ct^ (T, K) at time to of a call option with expiry T > to and strike 
K > is given by 



Ct,iT,K) = e-^*^(*^''*i?'[niax(5T - K,0)\Tto] 



(2) 



As argued below (see Section [TT^ under Assumption (|H]), the expectation in ^ 
is finite. 

Our main result is the following: 

Theorem 1 (Forward PIDE for call options). Let ipt be the exponential double 
tail of the compensator m(t, dy) 



Jl^ dx e^ floo '^") z <Q 
dx m{t, du) z > Q 



and define 



Xuy{z) ^ni^t{z)\St-^y] 



(3) 



(4) 



Under assumption ((Ht . the call option price {T,K) i— > Ct„{T,K), as a function 
of maturity and strike, is a solution (in the sense of distributions) of the partial 
integro- differential equation: 



dT ^ ' dK 2 dK^ JO 



+OC pi2 

y 



(T,dy)xT,Jln(- 



K 



y 

(5) 

on [to, cx3[x]0, oo[ with the initial condition: \/K > Cfo(to, K) = (Stg — K) + - 



Remark 1. The discounted asset price 

is the stochastic exponential of the martingale U defined by 
Ut^ f St dWt + I [{ey - l)M{dtdy). 



Under assumption (jH|, we have 
VT > 0, 

and \29l Theorem 9] implies that (St) is a f -martingale 



exp( ^(C/^[/=)T + (C/^C/'^)T 



< oo 
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The form of the integral term in ([S]) may seem different from the integral 
term appearing in backward PIDEs [24]. The following lemma gives an 
interpretation of XT,y{z) in terms of call payoffs, which casts this term into a 
more familiar form: 

Lemma 1. Let n{t, dz, y, lo) dt be a random measure on [0, T] x R x IR+ verifying 

/oo 
(e"^ A |zp)n(t,dz,y,w) < oo a.s. 

Xt,y{z) exponential double tail defined as 

Xt y{z)^l ^-f^'^'' (6) 

1 ir '^(*' y) ^ > ^ 

Then 

Jjiye^ - K)+ - e^{y - K)+ - K{e^ - l)l,>Kj]n(t, dz, y) = yxt,v (^In 
Proof. 

4[(ye- - K)+ -e^{y~ K)+ - K{e^ - l)l^>K}]n(i, dz, y) 
= 41(^6^ - ^)l^>in(f )} - - K)ly^K] - K{e^ - l)lj,>^}]n(i, dz,y) 
= /R[(ye^ ~ ^)lz>in(f )} + {K~ ye')ly>K}]nit,dz, y) 

• li K >y, then 

+ 00 

y(e-„e'"(f))n(t,dz,y) (7) 

• li K <y, then 
li<-<y} [(ye^^ - (|.)} + (A' - ye^)lyyK}]nit, dz, y) 

[{ye' ~ K) + {K- ye')]n(t, dz, y) + [K - ye']n{t, dz, y) 

ln(f) J-oo 

'"'"^(e'"(f^-e^Wt,dz,2/) 

-OO 

(8) 

Using integration by parts, xt,y can be equivalently expressed as 

M = / /-oo(e"-e")n(t,du,2/) z<0 
1 C(e"-e^)n(t,du,2/) z>0 
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Hence: 



Jjiye' - K)+ - e\y - K)+ - K{e/- - l)l,>K}]^(i, dz, y) = yxt,y (\n (^^)) 



□ 



1.2 Proof of main result 

In this section we present a proof of Theorem [TJ The proof is based on the 
Tanaka-Meyer formula [23l Theorem 9.43] and assumption pl| . 

Proof. We first note that, by replacing P by the conditional measure P|jro given 
J-Q, we may replace the conditional expectation in ^ by an expectation with 
respect to the marginal distribution p^{dy) of St under ^\j^t^- Thus, without 
loss of generality, we put io = in the sequel, consider the case where To only 
contains null sets and we denote Co(r, K) = C{T, K) for simplicity. 
([2]) can be expressed as 



C(r, K) = e- ^0 rW dt _ K)+ pS^dy) 



By differentiating with respect to K, we get: 



dC 
'dK 



poo 

Jk 



{T,dy)^e-fo^W^YT{dy) 



(9) 



(10) 



Let Lf = Lf{S) be the semimartingale local time of S at K under P (see (131 
Chapter 9] or [30l Ch. IV] for definitions). For ft, > 0, applying the Tanaka- 
Meyer formula to [St — K)^ between T and T + h, we have 

{St+h - K)+ = {St - K)+ + [^^'^ ls,_>K}dSt + \{L^^^ - L^) 

(11) 

+ {St-K)+ -iSt--K)+ -ls,_>K}^St. 

T<t<T+h 



As noted in RemarkfH the integrability condition (|H| implies that St — -^o ' '^'^St 
£{U)t is a martingale under P. So ([Ij can be expressed as dSt = r{t)St~dt + dSt 
where the discounted price : 

rT+h 



T+h rT+h 

^St->K}dSt = 



^St->K}dSt 



r{t)St-ls,_ 



IT JT JT 

where the first term is a martingale. Taking expectations, we get: 



>K} 



dt 



5/0^+" r{t) dtc(T + h,K)- efo dtC{T, K) 



-E 



= E 

'^T<t<T+ 



J^^\{t)St ls,_>K}dt+'^{L^^„ - Lf] 



J5t ~ K)+ - {St- - K)+ - \s,_>K}^St 
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Noting that St-lst_>K} = {St- — K)~^ + Klsf_yK}, we obtain 

pT+h ] pT+h 

/ rit)St-ls,_>K}dt = / rit)elo^(^)'^^ 
Jt Jt 

using Fubini's theorem and (1101) . As for the jump term, 

ET<t<T+hiSt - K)+ - [St^ - K)+ - ls._>K}A5t 



C(t,K)-K—{t,K) 



dt 



E 



= E 



/^+'' dt J m{t, dx) [St-e- ~ K)+ - {St- - K)+ - ls,_>K}St- (e" - 1) 



/^+'' dt J m{t, dx){{St-e=' - K)+ - {St- - K)+ 
{St- - K)+{e- - 1) - Kls,_yK}{e^- - 1)) 
/^+''diE[/ m{t,dx){{St-e^ - K)+ 
-e-{St--K)+ -Kls,_yK}{e^ -I) 



Applying Lemma [l] to the random measure m and its exponential double tail 



J m{t,dx){{St-e^-K)+-e%St--K)+~Kls,_>K}{e''-l)) ^ S 
holds, leading to: 



't- ^t,s 



E 



{St-K)+ ~{St--K)+ -ls,_>K}ASt 

T<t<T+h 
T+h 

dtE St-^tA- (In ( ^ 



T+h 



dtl 



i't,s,^ ( In ( ^ ) ) 1^*- 



T+h 



dtl 



St- 



St-Xt,St- [ In ( 



(12) 



Let (fi G C^(]0,oo[) be a test function. The occupation time formula (see [23l 
Theorem 9.46]) yields: 



T+h 



T+h 



dK^{K){Li^+^ -L^)^ I ^{St-)d[Srt = / dtcp{St-)SlSf 



{St) is a martingale, hence : E[S't] < oo. Since {Sr+h ~ K)^ < S'r+ft,, 
{St-K)+ < STAET<t<T+h{St-K)+-{St--K)+-ls,_>K}^St\ < 35* and 
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E 



St->K}dSt 



< oo then ([TT|) leads to : E [Ll^^/^ — L^] < oo; further- 
more, since ip is bounded and has compact support, one may take expectations 
on both sides and apply Fubini's theorem to obtain: 



/ dK^{K){L^^„-Lf^) = E f 
Jo J Jt 



pT+h 



^{St)St5tdt 



T 

T+h 

I 

T+h 



dm[p{St)s^j, 



E 



T+h 



dtipiSt-)S^_ait:St-) 



T 

oo rT+h 



(p{K)K^a{t, K fpf{dK) dt 

JT 

T+h poc fP'C 

dtef^^'^^"^ ^{K)K^a{t,Kf^{t,dK) 



where the last line is obtained by using (|T0|). Gathering together all the terms, 
we obtain: 



dK ip{K) po"'^" '^'C{T + h,K)- e^o '^^CiT, K) 

T+h i-co 

dtr{t)e^o^^^)d' / dK ^{K)[C{t,K) - K^{t,K)] 



T+h 



^KMt.Kf^{t,dK) 



fT+h 

+ / dKip{K) I dm 



St-xt,st- ( ( 



(13) 
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Dividing by h and taking the limit h —> yields: 



dKip{K)e 



J^r{t)dt 



poo 

Jo 



— iT,K)+r{T)C{T,K) 
fin 

r{T)C{T,K)-r{T)K—{T,K) 



dK(p{K)l 



K 

'S't-XT.ST- (In ( -J — )) 



dK^{K)e^o^W<it 



r{T)C{T,K)~r{T)K—{T,K) 



dK(p{K)e 



+ 00 c,2 



yi^'yT,dy)xT.,y ( In 



K 



JO dK^ ' - - \ \y 

Since this equality holds for any (p e C^(]0, cx)[, R), C(., .) is a solution of 



(14) 



dC 
df 



{T,K) 



2 

in the sense of distributions on [0, T] x]0, oo[. 



gK2iT,K)-riT)K—iT,K) 



□ 



2 Examples 

We now give various examples of pricing models for which Theorem [1] allows to 
retrieve or generalize previously known forms of forward pricing equations. 



2.1 Ito processes 

When iS* is an Ito process i.e. when the jump part is absent, the forward equation 
^ reduces to the Dupire equation [T7]. Theorem Theorem [T] gives a precise set 
of assumptions under which the Dupire equation holds, which do not require 
non-degeneracy conditions as in Gyongy [?]: 

Proposition 1 (Dupire PDE). Consider the price process S whose dynamics 
under the pricing measure P is given by: 



Define 



St = So+ I r{t)Stdt + / St5tdWt 
Jo Jo 



(15) 
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E 



exp 



6fdt 



< oo a.s. 



the call option price ^ is a solution (in the sense of distributions) of the partial 
differential equation: 



da 



to 



-r{T)K 



to 



(16) 



dT ' dK 2 

on [to, c»[x]0, c»[ with the initial condition: yK > C'to{to,K) — {St„ — iir)+. 

Proof. It is sufficient to take /i = in ([1]) then equivalently in ([S|). We leave the 
end of the proof to the reader. □ 



2.2 Markovian jump-difFusion models 

Another important particular case in the literature is the case of a Markov 
jump-diffusion driven by a Poisson random measure. Andersen and Andreasen 
derived a forward PIDE in the situation where the jumps are driven by a 
compound Poisson process with time-homogeneous Gaussian jumps. We will 
now show that Theorem [T] implies the PIDE derived in [2], given here in a more 
general context allowing for a time- and state-dependent Levy measure, as well 
as infinite number of jumps per unit time ("infinite jump activity"). 

Proposition 2 (Forward PIDE for jump diffusion model). Consider the price 
process S whose dynamics under the pricing measure P is given by: 



St — Sq- 



r{t)St-dt- 



StMt,St-)dBt 



+ 00 



St-iey-l)Nidtdy) 
(17) 

where Bt is a Brownian motion and N a Poisson random measure on [0, T] x 
R with compensator v{dz) dt, N the associated compensated random measure. 
Assume: 

cr(., .) is bounded {H'l) 

Then the call option price 

Ct,iT,K) = e-^*^(*'''*i?P[max(5T - K,0)\Tt,] 

is a solution (in the sense of distributions) of the partial integro- differential 
equation: 



dCto _^^T)K^Ct, 



dT 



2 dK^ 

Ct, (T, Ke-') - Ct, (T, K) - K{e-^ - 

on [to, oo[x]0, oo[ with the initial condition: > Ct„{tQ, K) 



dK 

v{dz) 



1) 



to 



dK 

{Sto 



(18) 



K) 
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Proof. We first note that, by replacing P by tlie conditional measure P|jro given 
J-Q, we may replace the conditional expectation in ([2]) by an expectation with 
respect to the marginal distribution p^{dy) of St under Thus, without 

loss of generality, we put to = in the sequel, consider the case where J-q only 
contains null sets and we denote Co(T, K) = C{T, K) for simplicity. 
By ^ differentiating in the sense of distributions with respect to K, we get: 



(19) 
(20) 



In this particular case, m(t, dz) dt = v{dz) dt and V't is simply defined by: 
Then (|4|) yields Xt,y{z) = xi^) — V-'(^)- Let now focus on the term 



dx vidu) z < 

dx v{du) z > 



+ 00 



,|^(T,dy)x(ln(^ 



in ([5]). Applying Lemma [T] 



e-/o"Kt)<i*pS(^y) y^[(ye- - K)+ - e\y - K)+ - K{e/' - l)ly^K}Hdz) 

poo 

/ e-fo^W^'pUdy)[{y-Ke-Y-{y-K)+~K{l-e-^)ly:,K^],.{dz) 



C{T, Ke-^) - C{T, K) - K{e-^ -I) — 



i>{dz) 



(21) 



This ends the proof. 



□ 



2.3 Pure jump processes 

We now consider price processes with no Brownian component. Assumption 
(IhI) then reduces to 



VT > 0, E 



exp 



dt / {e^ - lfm{tdy) 



< oo 



iH2) 
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and the forward equation for call option becomes 

It is convenient to use the change of variable: v — Iny, fc = InK. Define, 
c{k,T) = C(e'=,T). Then one can write this PIDE as: 

§ + "'^'i = /r (§ - 1) - 

where xt,v is defined by: 
with: 

In the case, considered in [S], where the Levy density my has a deterministic 
separable form: 

my (t, dz, y) dt — a{y, t) k{dz) dz dt (24) 
Equation (|23p allows us to recovei0 equation (14) in 9 : becomes 



where k is defined as the exponential double tail of k{u) du, i.e: 
k{z) = 



dx k{u) du z < 

J^°° dx k{u)du z > 



The right hand side can be written as a convolution of distributions: 

|£ + „r)|^[a.(0 (g-l)].. wh™ (25) 

g{u)^e-^"K{u) aT{u)^a(e^'',T) (26) 

Therefore, it implies that from the knowledge of c(., .) and a choice for k(.) 
we can recover ax hence a (.,.). As noted by Carr et al. |9j, this equation is 
analogous to the Dupire formula for diffusions: it enables to "invert" the struc- 
ture of the jumps-represented by a- from the cross-section of option prices. 
Note that, like the Dupire formula, this inversion involves a double deconvolu- 
tion/differentiation of c which illustrates the ill-posedness of the inverse problem. 



^Note however that the equation given in [5] does not seem to be correct: it involves the 
double tail of k(z) dz instead of the exponential double tail. 
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2.4 Time changed Levy processes 

Time changed Levy processes were proposed in [5] in the context of option 
pricing. Consider the price process S whose dynamics under the pricing measure 
P is given by: 

St - exp (Xt + / r{t) dt) Xt = exp (LtJ Tt = j 9,ds (27) 
Jo Jo 

where Lt is a Levy process with characteristic triplet (6, N its jump 

measure and {9t) is a locally bounded positive semimartingale. We assume L 
and 9 are J^t-adapted. 

Xt = (e~ -I'd St) is a martingale under the pricing measure P if exp (Lt) is 

a martingale which requires the following condition on the characteristic triplet 
of (Lt): 

b+la^+ [ {e'-l-z l|,|<i)j^(dy) = (28) 

Define the value Ct^, (T, K) at time to of the call option with expiry T > to and 
strike if > of the stock price (St): 

Ct„{T,K) - e-fo-WE^[ma^{ST - K,0)\gt„] (29) 
Proposition 3. Define 

a{t,x)^E[9t\Xt-^x] 
and X the exponential double tail of v{du) 



^ dx v{du) z < 

dx v{du) z > 



Then under assumption {H2), the call option price Cto ■ {T,K) H- Cta{T,K) 
at date tg, as a function of maturity and strike, is a solution (in the sense of 
distributions) of the partial integro- differential equation: 

dT ra^-^^^^)^ Qj^+ 2 dK^ ^^^^ 

on [i, oo[x]0, oo[ with the initial condition: Vif > Cta{to, K) = {Stg — K)+- 
Proof. Using Aj Theorem 4] , {Lt^ ) writes 

Lt, ^Lo + /o cT^/TsdEs + /o b9sds 

+ /o ^« /|.|<i ^^^ds dz) + 9, 4,^^ zN{ds dz) 
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where N is an integer-valued random measure with compensator v{dz) dt, N its 
compensated random measure. 
Applying the Ito formula yields 



Xt 



+ ^ Xs-dLr^ + i / ds 



"^Xs-Xs--Xs-AL 
1 , 



= Xo 



Xs- 







ds 



Xs-aJBs dBs 



+ / Xs-Os / zN{dsdz)+ / Xs-0s / zN{dsdz) 



Xs-Osie^ - 1 - z)N{dsdz) 



Under assumption (i?2a)i /(^^ — 1 — z l|2|<i)i^((iz) < oo, hence: 



Xt 



1 



Xo+ / Xs- 

rt 

' Xs-Osie^ - l)N{dsdz) 



be, + + I {e'-l- zli,i<,)dsiy{dz) 



ds+ I Xs-aVedB, 
'o 



"'R 



= Xn 



Xs-cr^/e'sdBs+ / / Xs-0s{e"- -l)N{dsdz) 



and (St) may be expressed as: 

St = So+ f Ss-r{s)ds+ f Ss-(Jy%dBs+ f f Ss-Os{e/ - l)N{ds dz) 
Jo Jo Jo Jm 

(St) is now in the suitable form ([T]) to apply Theorem [TJ which yields the 
result. □ 



2.5 Index options in a multivariate jump-diffusion model 

Consider a muhivariate model with d assets: 

5-^ = 5*^+ [ r{t)Sl-dt+ [ St-SldWi+ [ [ Sl-ie^' -l)N{dtdy) 
Jo Jo Jo Jr'' 

where is an adapted process taking values in M representing the volatility of 
asset i, W is a. d-dimensional Wiener process, A'' is a Poisson random measure 
on [0,T] X M'* with compensator v{dy)dt, N denotes its compensated random 
measure. 

The Wiener processes W"^ are correlated: for all 1 < < d, {W'',W^)t = 
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pijt, with pij > and pu — 1. 

An index is defined as a weighted sum of the asset prices: 

d 

The value Ct^ (T, if) at time of an index caU option with expiry T > to and 
strike K > is given by 



Ct,{T,K) = e-^*o'-(*)'^*£;P[max(JT - if,0)|J-tJ 



(32) 



The foUowing result is a generalization the forward PIDE studied by Avellaneda 
et al. [3J for the diffusion case: 



Theorem 2. Forward PIDE for index call options. Assume 



cxp(i/o^||<5,|prft 



< oo 



VT > E 

/r<j(1 A II2/II) i'(dy) < oo a.s. {A2b) 

Let k{.,t,dy) be the random measure: 

'Y.l<^<a-l^^Sl-ey^+waSteV' 



(33) 



k{t,dy) 



In 



i-{0} 



v{dyi,--- ,dyd-i,dy) 



and rjtiz) its exponential double tail 
Vtiz) = 

Define: 



J^°° dx k{t, du) z > 



(34) 



a{t,z) 



\ 



E 



E [,74 (z) |/,_ = y] 



(35) 



(36) 



T/ie index call price (T, K) 1— >■ Ct„ (T, if) , as a function of maturity and strike, 
is a solution (in the sense of distributions) of the partial integro- differential 
equation: 



dC_ 
df 



-r{T)K 



dC a{T,KY d^C 



dK 



y^iT,dy)xT,y [ In 



(37) 

on [io, oo[x]0, oo[ with the initial condition: Vif > Cto{to,K) = (it^ — if)+. 
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Proof. {Bt)t>o defined by 



1/2 



is a continuous local martingale with quadratic variation t: by Levy's theorem, 
B is a Brownian motion. Hence / may be decomposed as 



^=l U,j=l 



(38) 



r [ y^w,SlJey'-l)Nidtdy) 



The essential part of the proof consists in rewriting (It) in the suitable form ([T]) 
to apply Theorem [T] 
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Applying the Ito formula to In {It) yields: 

In {It) = ln{Io) + J^ 7~ + ^ / ^ ( |" ^" 
+ Yl ln(/t-+A/t)-ln(/t-)--^A/t 

= In (7o) + £ ^r{t) - 2^ E ^^^^'^''^ ^t^i Si-Sf^j j 

Jo Jk'' \ ^t- I 



•t-SL 1 at 



dt 



In (Jo) + r (r(t) - ^ I E ""^"^JPij ^tH SlSl I 



+ r f lnl ^'-'-'''''''^''''' ]N{dtdy) (39) 



JM<i 



The last equality is obtained since 

ity property of the logarithm ( one recalls that J2i<i<d it ~ 
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Holder inequality: 



In 



Ei<i<dW'i^l-e* 



< 



EWi, 
~li 



SI 



(40) 



l<i<d 



1/2 



(42) 



l<i<d 



, l<i<<i 



hence the functions y ^ In 



WiSl_eVi 



and y 



5Di<i<<i ■WiSl_e^ 



are in- 



tcgrabk; with respect to ^{dy) by assumptions (^26) and (^3;,). We furthermore 
observe that 



In 



y{dy) < 00 a.s. 



2 In 



/^ £i<.<d"'.gj-^''- M 



(43) 



J||?/||>i 



(dj/) dt < 00 a.s. 



Similarly, since (^25) and (^sb), for all 1 < i < d, J (e^' — 1 — l^y.^^ciVi ) '^(c^y) < 
00 and In (5^) rewrites: 



In (5^) = In (5§) + {r{t) - i(5j)' - iZo i^'' ' 1 - l|..l<i2/i ) K^^) '^i) 



dt 



Define the d-dimcnsional martingale Wt = (W/,--- ,W/ ^,-Bt). For all 1 < 
i,i<d-l, {W\W^)t = pijt; 



{W\B)t 



Ell WjPijSl_5l 



1/2 



Define 



pii 



172 



Pl,d-1 
Pd-l,d-l 



£j=l"'jPd-l,,.9;_ 

1 



There exists a standard Brownian motion (Zt) such that Wt = AZt where A is 
a d X d matrix verifying 
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Define Xt = (in (S*-^), • • • , In (5'^"^), In {It)) ; 



r{t)-USl)'-JZ^{ey^-l-y^) v{dy) 



r{t)-kistr~.r^i^''- 



1 - Vd-i) iy{dy) 



and 



( y\ 

My) = 



1 - In 



El 



In 



Vd-i 



\ 

then (Xt) may be expressed as: 

Xt = Xo+ I Ptdt+ [ dtAdZt 
Jo Jo 

Let us now compute the inverse ^ of tp: 



JR<* 



My)N{dtdy) (44) 



(i,w, (?/i,--- ,yd)) {4>t{i^,y)i,--- ,<i)t{i^,y)d) 

For all t in [0, r],for all y 

(l)t{uj,ipt{uJ,y)) = tjjt{i^,(pt{i^,y)) = y. 
Obviously, for 1 < i < (d — 1), {4>t{y))i = yi- For i = d: 



(45) 
(46) 



MMy))d = In 



= In 



jl<i<d- 



Hence: 



c(?/) = \yi,--- ,J/d-i,ln 



(47) 
= yd (48) 



WdSt 
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Observe that ipt{-,0) = 0, </> is predictable, and is difFerentiable with 

Jacobian matrix \/y(j)t{y): 



( 







e''<J/t--El<i<d-l™.Sj'_e« e!'<i/t_-El<i<d-l"'iS;_e«i / 



(-0, v) satisfies [H Assumption (H3)]: since (^2b), for all T and for all t in [0, T], 
(lA||^t(.,y)|P)Kd2/)dt 



E 



JM.!^ 
T 



< 



1 A I + • • • + yj_, + In 
1 A {2\\yf)uidy)dt < oo 



dt 



/O JR'' 

Define j/^, the image of v by 0: 

e S(K'^ - {0}) i, A) = iy{M^, A)) 

Applying [H Lemma 2], Xt may be expressed as: 

XT = Xa+ [ I3tdt+ [ 5tdZt+ I [ yM{dtdy) 



(49) 



where M is an integer-valued random measure (resp. M its compensated ran- 
dom measure) with compensator 

//(cj; dt dy) — m{t, dy; w) dt 

defined via its density with respect to v^: 



dfi 



1- E ^^r-^^" 



l<i<d-l 



(50) 



Considering now the d-th component of Xt, one obtains the semimartingale 
decomposition of In (It): 



ln(/T) - ln(/o 



mwjp.j 5\5l Sl_Sl_ 

T,l<^<dW^Sl_ey' 



- 1 - In 



I'idy) ] dt 



T 









dBt + 
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where K is an integer- valued random measure on [0, T] x R with compensator 
k{t, dy) dt defined by: VB € S(R - {0}), 



k{t,B) = ^i{t,W^-^ X B) 
and K its compensated random measure. Therefore 



(51) 



r(t)/t- dt- 



The normahzed volatiUty of It satisfies for all t in [0,r], 



I?- 



and 



Hence: 



In 



< 



1 



2 7o 

1 /-^ E'.^i^.^.p^.'^l'^^gL^^- 



Js.'' 



- l)It-K{dtdy) 



- lfk{t, dy)dt 



2 7o 



It 



dt 



Ei<i 



<(i-l ■"'^'-'t 



siey^ + wdStey 



JM'' 
d 



-1 v{dyi,--- ,dyd-i,dy)dt 



< 



Using assumptions (Au), {A2b) and (A^h), the last inequality implies that (It) 
satisfies (|H]). Hence (It) is now in a suitable form to apply Theorem [U which 
yields the result. We leave the end of the proof to the reader. □ 



2.6 Forward equations for CDO pricing 

Portfolio credit derivatives such as CDOs or index default swaps are derivatives 
whose payoff depends on the total loss Lt due to defaults in a reference portfolio 
of obligors. Reduced-form top-down models of portfolio default risk [5D1 HU [211 
nil I32j represent the default losses of a portfolio as a marked point process 
{Lt)t>o where the jump times represents credit events in the portfolio and the 
jump sizes ALt represent the portfolio loss upon a default event. Marked point 
processes with random intensities are increasingly used as ingredients in such 
models [2Ql EH [Ml ISH |32] . 
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In all such models the loss process (represented as a fraction of the portfolio 
notional) may be represented as 



Lt= [ [ xM {ds dx) 
Jo Jo 



where M{dt dx) is an integer- valued random measure whose compensator fi{dt dx; uj) 
m{t, dx; to) dt has finite mass 

Xt{uj) = / m{t,dx;uj) 
Jo 

Nt represents the number of defaults and A represents the default intensity i.e. 
the (random) jump intensity of the point process Nt = M{[0,t] x [0, 1]). Denote 
by Ti < T2 < .. the jump times of N. The marked point process L may also be 
represented as 

Nt 

Lt = ^ Zk 

k=l 

where the "mark" Zk taking values in [0, 1] is distributed according to 

mx{t, dx; w) 



Ft{dx]uj) 



At(w) 



Note that the percentage loss Lt belongs to [0, 1], so AX^ g [0, 1 — For the 
equity tranche [0, K], we define the expected tranche notional at maturity T as 

Ct,{T,K)^n{K~LT)+\Fto\ (52) 

As noted in [11] , the prices of portfolio credit derivatives such as CDO tranches 
only depend on the loss process through the expected tranche notionals. There- 
fore, if one is able to compute Ct^ {T, K) then one is able to compute the values 
of all CDO tranches at date Iq. In the case of a loss process with constant loss 
increment, Cont and Savescu [121 derived a forward equation for the expected 
tranche notional. The following result generalizes the forward equation derived 
by Cont and Savescu [12] to a more general setting which allows for random, 
dependent loss sizes and possible dependence between the loss given default and 
the default intensity: 

Proposition 4 (Forward equation for expected tranche notionals). Define the 
integer-valued random measure Myidtdy) with compensator mY{t,dy,Yt-) dt 
defined by : 

VAeS([0,l]), mYit,A,z) = E[mxit,A,.)\Lt- = z] (53) 
and the effective default intensity 

A^(t,z)- / mY{t,dy,z) (54) 
Jo 
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The expected tranche notional {T,K) h> Ctf,{T,K), as a function of maturity 
and strike, is a solution of the partial integro- differential equation: 



dC 



dT Jo dK^ 



K Q2 



{T,dy) 



K-y 



(K-y-z) my (T, dz, y) - (K ~ y)X'' {T, y) 



(55) 

on [to, oo[x]0, 1[ with the initial condition: yK £ [0,1] Cto{to,K) ~ {Lfg — 
K)+. 

Proof. By replacing P by the conditional measure given J-q, we may re- 
place the conditional expectation in (I52p by an expectation with respect to the 
marginal distribution pxidy) of Lt under Pij^jq- Thus, without loss of general- 
ity, we put to = in the sequel and consider the case where J-q only contains null 
sets. Observe that the definition of my consists in defining the compensator of 
the mimicking process of Ly- 

Define the operator Ht as follows: for / in C^{[0, 1]), 

HTf{y)= I \f{y + z)~f{y)]my{T,dz,y) 



Applying [U Theorem 4] leads to the Kolmogorov forward equation for ((if ), P|jr„) 

— - Hri..pT, 

where is the L^-adjoint of Ht- Hence for y — > {K — y)+ 



Hr.iK-y) 

UK>y, 
HT.{K-y)- 



+ _ 



[{K-y- z)+ -{K- y)+]my (T, dz, y) 



[{K -y - z)lK-y>z} -{K - y)lK-y>o}]'myiT, dz, y) 



\{K-y-z)-{K-y)]my{T,dz,y)+ [ \{0 - {K - y)]my{T,dz,y) 

JK^y 



K-y 



zmy{T, dz,y) - {K-y) 



my{T, dz,y) 



K-y 



my{T, dz,y) 



K-y 



{K-y-z) my{T,dz,y)-{K-y)X^{T,y) 
For K < y, Ht-{K — y)+ = 0. ([5^ can be expressed as 

C{T,K)= [\k - y)+ pT{dy) (56) 
Jo 
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By differentiating with respect to K, we get 
dC n , , . d^C 



dK 



Pridy) 



K 



{T,dy) ^pT{dy) 



dC_ 

dT 



{K - y)- 



dpT 
dT 



(dy) 



f {K-y)+ H^.pTidy)^ f pridy) Ht . {K - y)+ 
Jo Jo 



K 



Pridy) 



K-y 



(K-y-z) mviT, dz, y) ~ [K ~ y)X^iT, y) 



Then 

dC_ _ 
dT ~ 



K ^2 



iT,dy) 



K-y 



(K-y^z) mviT, dz, y) ~ {K - y)X'' {T, y) 



□ 



In [12], loss given default (i.e. the jump size of L) is assumed constant 
5 = {1 — R)/n: the marks are then deterministic and equal to 6 : Lt — SNt 
and one can compute C(T, K) using the law of Nt- Setting io = and assuming 
as above that J-t^ onl 

C{T, K) = K[{K - Lt)+] = E[(fc 6 - Lt)+] = (5IE[(fc - Nt)+] = S Ck{T) (57) 

The compensator of Lt is Xt es{dz) dt, where es{dz) is the point mass at 
the point S. The effective compensator becomes: mY(t,dz,y) = E[Xt\Lt^ = 
y]es{dz)dt = X^ {t,y) es{dz) and the effective default intensity is X^{t,y) = 
E[Xt,\Lt^ = y]. 

If we set y = j6 then : X^{t,jS) = E[Xt\Lt- = jS] = E[Xt\Nt- = j] = aj{t) 
and pt{dy) — X]^=o 1ji^)^js{dy) with the notations in il2 . 
Let us focus on (l55p in this case. We recall from the proof of Proposition|3]that: 



dC 
df 



(T, kS) 



pridy) HT.{k5-y)+ 







'pridy) f \ik5~y-z)+ -ik5-y)+]X''iT,y)esidz) 
Jo 
1 







pridy) X^ (T, y) [(M - y - 5)+ - (M - y)+] l,-<i_,} 



-^^(IjiT)ajiT) lj<fc-i 

3=0 



(58) 
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This expression can be simplified as in [T^l Proposition 2] , leading to the forward 
equation 

dC (T) k—2 

= afe(m_i(r)-afc_i(r)Cfe(T)-^Q(T)[a,+i(T)-2a,(T) + a,_i 

fc-2 

= [ak{T) - ak-i{T)]Ck-i{T) - ^(V2a),Q(T) - ak-i{T)[Ck{T) - Cu- 
Hence we recover [121 Proposition 2] as a special case of Proposition HI 
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